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Noise  of  Mode-Locked  Lasers 

H.  A.  Haus,  Life  Fellow.  IEEE,  and  Antonio  Mecozzi 


Abstract— A  theory  of  noise  in  mode-locked  lasers  is  devel¬ 
oped  that  applies  to  additive  pulse  mode-locked  and  Kerr  lens 
mode-locked  systems.  Equations  of  motion  are  derived  for  pulse 
energy,  carrier  linewidth,  frequency  pulling,  and  timing  jitter. 
The  effect  of  gain  fluctuations,  mirror  vibrations,  and  index 
fluctuations  are  determined.  Measurements  that  can  determine 
all  four  fluctuation  spectra  are  described.  Experimental  data  in 
the  literature  are  compared  with  theory. 

I.  Introduction 

IN  recent  years,  many  solid-state  lasers  have  been  mode 
locked  successfully,  using  the  additive  pulse  mode 
locking  (APM)  principle  [l]-[6].  It  has  been  shown  that, 
under  appropriate  approximations,  the  theory  of  APM  is 
congruent  with  that  of  saturable  absorber  mode  locking, 
supplemented  by  self-phase  modulation  (SPM)  and  group 
velocity  dispersion  iGVD)  17].  Another  system  of  passive 
mode  locking  has  been  dubbed  Kerr  lens  mode  locking 
[8]-[10]  (KLM).  It  utilizes  self-focusing  to  produce  an 
intensity  dependent  loss.  It  is  described  by  the  same  the¬ 
ory. 

Some  of  the  systems  are  already  incorporated  in  com¬ 
mercial  products.  As  time  proceeds,  different  schemes  will 
be  selected  according  to  their  ability  to  produce  stable 
short  pulses  with  the  required  duration  and  intensity.  It  is 
likely  that  the  susceptibility  to  environmental  noise  will 
differ  among  the  different  systems.  One  purpose  of  this 
paper  is  to  predict  the  noise  characteristics  as  governed 
by  specified  pump,  mirror-position,  and  beam  axis  fluc¬ 
tuations.  Another  purpose  is  more  fundamental.  To  the 
authors’  knowledge,  an  inclusive  theory  of  noise  in  mode- 
locked  systems  with  saturable  absorber  action,  SPM  and 
GVD  does  not  exist.  The  reason  for  this  is  that  perturba¬ 
tion  theories  developed  for  Hamiltonian  systems  do  not 
apply;  the  system  is  not  self-adjoint.  In  this  paper  we  de¬ 
velop  an  analytic  theory  of  adiabatic  pulse  evolution  in 
the  presence  of  noise.  We  take  advantage  of  the  fact  that 
the  shortest  pulses  are  achieved  when  they  are  transform 
limited  and  secant  hyperbolic  in  temporal  shape.  If  GVD 
and  SPM  are  important,  they  are  perturbed  solitons.  In 
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this  limit,  a  perturbation  theory  is  feasible,  closely  related 
to  that  of  soliton  perturbation  theory  [II].  Yet  even  in  the 
case  when  the  simple  saturable  absorber  mode-locking 
equations  apply  [12],  the  perturbation  equations  can  be 
used,  qualitatively,  with  proper  interpretation  of  the  pa¬ 
rameters. 

In  Section  II  we  review  the  master  equation  of  saturable 
absorber  mode  locking  in  the  presence  of  (negative)  GVD 
and  summarize  its  solution.  Section  III  investigates  a  lin¬ 
earized  equation  for  the  perturbation  from  the  steady-state 
solution  as  driven  by,  as  yet,  unspecified  noise  sources. 
The  perturbed  pulse  is  described  by  four  parameters,  the 
amplitude  perturbation,  and  the  perturbations  of  phase, 
frequency,  and  timing.  In  addition,  there  is  a  continuum 
that  is  not  part  of  the  perturbed  pulse.  In  order  to  project 
out  the  amplitudes  of  the  perturbation  from  a  specified 
complex  amplitude  function,  one  uses  solutions  of  the  ad¬ 
joint  system  which  are  orthogonal  to  the  solutions  of  the 
linearized  equation  of  motion  (Section  IV).  Using  the  or¬ 
thogonality  properties,  equations  of  motion  are  written  for 
the  four  pulse  perturbation  parameters.  These  are  linear 
coupled  first  order  differential  equations  in  time,  driven 
by  noise  sources. 

The  analysis  ignores  the  coupling  to  the  continuum. 
This  is  legitimate  as  long  as  the  continuum  does  not  react 
back  appreciably.  One  may  develop  plausibility  argu¬ 
ments  that  the  approximation  is  valid.  Section  V  identifies 
the  noise  sources  starting  from  different  realizations  of 
APM  systems  and  KLM  systems  and  expresses  the  noise 
sources  in  terms  of  physical  perturbation  spectra.  Section 
VI  evaluates  the  theoretical  spectra  and  the  correlation 
functions  of  the  four  observables.  Section  VII  determines 
ways  of  measuring  the  spectra  of  each  of  the  four  pulse 
parameters.  Section  VIII  compares  experimental  data  with 
theory. 


II.  The  Master  Equation  and  its  Solution 
The  master  equation  for  the  nonlinear  pulse  evolution 
as  a  function  of  the  slow  time  variable  T  (expressing  the 
evolution  over  many  cavity  roundtrips)  is  [7] 


-l  +  g 


±ii\ 

(1,  dr  +  (l2gdt2/ 


+  (7  —  ;*)|a|2 


a  +  TRS(r,  T). 


(-•; 


Here,  TR  is  the  round-trip  time;  a(T,  t)  is  the  (electric 
field)  amplitude  in  the  resonator,  a  function  of  the  short¬ 
term  time  variable  t  as  well  as  of  T,  a  time  variable  on  the 
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scale  of  many  cavity  round-trip  times;  /  is  the  incremental 
loss,  g  the  incremental  gain  upon  one  pass,  is  the  gain 
bandwidth;  y  is  the  effective  saturable  absorber  action  that 
may  be  produced  by  APM  action,  6  is  the  Kerr  phase 
modulation  coefficient,  D  is  the  group  velocity  dispersion 
coefficient;  5(r,  T)  is  the  noise  source  [10).  The  operator 
multiplying  the  gain  g  expresses  the  dependence  of  the 
group  velocity  upon  the  gain  and  the  “diffusion”  of  a 
pulse  passing  the  gain  medium  because  of  the  finite  gain 
bandwidth. 

The  gain  is  assumed  to  have  a  long  relaxation  time,  so 
that  the  change  of  gain  upon  one  pass  can  be  ignored.  The 
gain  is  saturated  by  many  pulses  in  succession.  It  is  re¬ 
lated  to  the  small-signal  gain  g„  by 


where  Ps  is  the  saturation  power.  The  saturation  law  (2) 
neglects  dynamics  related  to  gain  relaxation.  Thus  effects 
that  depend  critically  on  the  gain  relaxation  time  are  not 
considered.  The  amplitude  a(T,  t)  is  so  normalized  that 
|a|:  is  equal  to  the  power.  The  solution  of  (1)  is  [12],  [13] 


a(t,  T)  =  A0  sech 


1 

-(/  -  O 
r 


ill  * jS) 


exp 


(3) 


where 


j  +  g  -  /  + 


(i  +  n 3) 


2  r 


8  ,  -r, 

M  J? 


=  o  (4) 


l*+jD 


(2  +  3/0  -  /32)  =  (7  -  jb)A\  (5) 


and 


t„  = 


n  J* 


T. 


(6) 


Here,  (3  is  the  chirp  parameter;  0  is  the  phase  of  the  pulse, 
an  arbitrary  parameter  in  the  noise-free  solution,  yet  a 
driven  quantity  when  noise  is  introduced;  ^  is  the  phase 
shift  produced  per  pass,  and  determined  from  the  solution 
of  (4).  Upon  multiple  passes,  this  phase  shift  per  pass 
translates  into  a  net  frequency  shift. 

Equation  (5)  shows  that  the  area  A0t  of  the  absolute 
magnitude  of  the  pulse  amplitude  is  fixed  for  a  fixed  set 
of  parameters  of  the  system.  The  perturbation  analysis  can 
proceed  from  two  starting  points.  One  is  the  original 
equation  of  fast  saturable  absorber  mode  locking,  in  the 
absence  of  GVD  and  SPM  [12],  [14]  with  5  =  D  =  0. 
One  may  find  a  set  of  eigenfunctions  of  the  inphase  and 
quadrature  perturbations  and  expand  the  noise  sources  in 
this  orthogonal  set  of  eigenfunctions.  Each  of  the  ampli¬ 
tudes  of  the  expansion  obeys  a  first-order  differential 
equation  in  T.  SPM  and  GVD  can  then  be  treated  as  ad¬ 
ditional  pefurbalions.  Wc  shall  not  pursue  this  approach 


here,  because  systems  with  genuine  “physical”  fast  sat¬ 
urable  absorbers  are  not  of  practical  importance,  because 
of  the  difficulty  in  finding  suitable  absorbing  materials. 
Of  much  greater  practical  importance  have  become  sys¬ 
tems  that  construct  artificial  saturable  absorbers  by  APM 
or  KLM  action.  Such  systems,  generally,  have  apprecia¬ 
ble  SPM  action  compensated  by  negative  GVD.  For  this 
reason,  the  pulses  are  soliton-like,  yet  not  solitons;  they 
are  solitary  pulses.  Contrary  to  regular  soliton  perturba¬ 
tions,  amplitude  and  frequency  perturbations  of  mode- 
locked  solitary  pulses  do  not  persist,  but  “shed”  some  of 
their  energy  into  the  continuum.  This  is  a  consequence  of 
the  gain  saturation  and  gain-bandwidth  limitation  not  en¬ 
countered  with  solitons  (nonlinear  Schrodinger  equation. 
NLSE).  The  feedback  of  the  continuum  radiation  back 
into  the  pulse  is  of  second  order  in  this  coupling  to  the 
continuum  and  will  be  neglected.  This  approximation  will 
be  confirmed  by  computer  simulation  later  on  in  the  pa¬ 
per. 

The  chirp  parameter  (j  is  zero  when  the  gain  dispersion 
and  APM  action  are  related  to  the  SPM  and  GVD  by 


g/ti]  =  7 
-D  6 


n- 


(7) 


Since  6  >  0  for  the  usual  Kerr  effect  and  y  >  0  for  sta¬ 
bility,  this  is  the  case  of  negative  GVD.  This  adjustment 
leads  to  transform  limited  pulses.  Henceforth  we  shall 
concentrate  on  this  special  case.  One  has 


0 


a;t- 


(8) 

(9) 


8  ~  I 

Substitution  of  these 
a„(t,  T)  =  A„  sech 


U>|  _ 

2 

--A2 

(10) 

7 

2 

i  g 
r2n; 

=  -- A 2 

1A° 

(ID 

relations  into  (3)  gives 


(12) 


T 


i  lm 

A0\  6 


(13) 


III.  The  Linearized  Perturbation  Equation 
Noise  is,  of  course,  handled  by  perturbation  theory.  If 
the  noise  is  small  compared  with  the  signal,  the  equations 
can  be  linearized.  In  the  adiabatic  case,  the  pulse  param¬ 
eters  are  changed  by  the  noise,  resulting  in  new  pulses 
that  obey  the  steady  state  equations.  The  area  of  the  pulse 
(amplitude-width  product)  stays  invariant.  An  adiabatic 
perturbation  of  the  pulse  changes  its  amplitude,  phase, 
frequency  and  timing. 
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We  make  the  ansatz: 

a(t,  T)  =  [a,(t  -  /„)  +  Aa(t  -  t0,  T)] 


(  5  ,2  T 

i*°T' 

where 


as(t)  =  A0  sech  (15) 

Substitution  of  the  expansion  (14)  into  (1)  and  disregard 
of  terms  higher  than  first-order  in  A  a  leads  to: 

t'Ttmu  ~ T) 

=  (A  +  |0| 

+  2taUl  ~  O ]  +  (M  -  O&0* 

+  TRS(t,  T).  (16) 

We  suppressed  the  indication  of  the  independent  variable 
t  -  t„.  The  parameter  gs  has  been  defined  as: 

_  So _ 2tA„ 

oj  n  T  /-  7  •  \  / 

{■+^i 

It  is  convenient  to  introduce  the  pulse  energy  w„  as  a 
new  parameter: 

»o=\  dt  a]{t)  =  2A2oT  =>A0  =  (18) 

We  write  (16)  in  the  formal  way: 

Tr  ~  Aa(t  -  t„)  =  A,(Aa)  +  TRS(t-  ta,  T)  (19) 

where  the  operator  A,  has  been  defined: 

6 

A,(An)  =  (fi  -  j)  -  -Al  +  |D |  —2 


IV.  Approximate  Solution  of  Perturbation 
Equation 

Small  changes  of  pulse  energy,  phase,  timing,  and  fre¬ 
quency  (—p  —  u>  ~  (ji0,  the  deviation  from  the  angular 
carrier  frequency  (11]  are  derivatives  of  the  steady  state 
solution  (12): 

fw(t)  =  =  —  1  ~  -  tanh  (-)  as(t)  (21) 

avv’  Wn  7  \  7 ) 


feU)  =  ~  =  jas(t) 

.  dau  1 
fit)  =  -  ~  tanh 

Otn  7 


+  2ba]  0  ~  K)  &a(t  ~  t„) 

+  (p  -  j)Sa?(t  ~  t„)Aa*(t  -  t„) 

~  ^  “-u  ~  ° 

•  j  dt  a,(Aa  1  A a*).  (20) 

In  the  next  section,  we  express  Aa  in  terms  of  the  pertur¬ 
bations  of  the  four  pulse  parameters. 


ip  =  i  --  ta,(t). 
w. 


We  expand  the  perturbation  in  terms  of  these  four  func¬ 
tions.  The  remainder  is  the  continuum  that  is  not  part  of 
the  pulse. 

A a  =  fK(t  -  t„)Aw(T)  +f,{t  -  to)A0(T) 

+  fP(t  ~  t„)Ap(T) 

+  f(t  -  ta)At(T)  +  A ac(t  -  t„,  T).  (25) 

Soliton  perturbation  theory  starts  with  the  above  expan¬ 
sion  and  follows  the  four  perturbation  amplitudes  Aw,  A 0, 
A p,  and  At  as  they  evolve  in  Til  1]  according  to  the  non¬ 
linear  Schrodinger  equation  (NLSE).  In  this  case.  Aw  and 
Ap  are  invariants,  AO  changes  due  to  the  Kerr  phase  shift 
change  caused  by  Aw,  and  At  changes  because  the  group 
velocity  changes  with  a  change  of  frequency  (Ap  ).  The 
master  equation  (1)  contains  two  additional  effects, 
namely  gain  dispersion  (GD)  and  APM  action.  As  a  con¬ 
sequence  of  these  effects.  Aw  and  Ap  do  not  remain  in¬ 
variant.  If  these  effects  are  small,  i.e.,  if  the  parameters 
p  defined  in  (7)  is  much  less  then  unity,  GD  and  APM 
can  be  treated  as  perturbations  on  the  NLSE.  The  cou¬ 
pling  to  the  continuum  of  the  changes  of  Aw  and  Ap  is  of 
order  p,  the  coupling  of  the  continuum  back  to  the  pulse 
is  of  order  p  square.  Hence  one  may  use  soliton  pertur¬ 
bation  theory  to  evaluate  the  four  pulse  perturbation  pa¬ 
rameters  to  first  order  in  p. 

Soliton  perturbation  theory  associates  a  set  of  orthog¬ 
onal  adjoint  functions  with  the  four  functions  and  the  con¬ 
tinuum  (11).  We  choose  them  so  that  the  products  inte¬ 
grate  to  unity,  unless  they  are  orthogonal  and  the  integrals 
are  zero. 


jfw(t)  =  2a,  (0 

(26) 

/•(/)  =  2/—  1 

-  tanh  ^  j  ajt) 

(27) 

fit)  =  —  ta,(t) 

(28) 
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fp(t )  =  j  [—  tanh  -  )  as(t).  (29) 

\woT  Tj 

One  may  use  the  orthogonality  property: 

Re  j  dtff(fjfj(t)  =  bu  (30) 

to  project  out  the  coefficients  of  the  perturbation  expan¬ 
sion: 


Aw(T)  =  j 

j  dt  [ft(t)Aa(t)  +  fw(t)Aa*(t)] 

(31) 

AW  =  j  | 

[  dt  [f*(t)Aa  (/)  +  fe(t)Aa*(t)} 

(32) 

A p(T)  =  i  1 

« 

j  dt  [fZ(t)Aa(t)  +  fp(t)Aa*(t)) 

(33) 

A  t(T)  =  ^  1 

f  dt  [fZ(t)Aa(t)  +  f,(t)Aa*(t)]. 

(34) 

With  the  aid  of  these  projection  operations,  one  may  de¬ 
rive  equations  of  motion  of  the  pulse  perturbation  param¬ 
eters. 

By  making  use  of  the  identity 


5  2  „  d2  , 

-A„  +  D~2  +  8as 


as  =  0 


(35) 


the  following  identities  can  be  proven  for  the  operator  A, 
defined  in  (20) 


AAfAD] 


{fi-j)6A20  -  gs 


( 


1  -  — - 
Qgdt 


A,[M‘)  1  =  0 

A,[fp(t)]  =  2\D  |  (1  +jn)f 
A,[M  01  =  0. 


as(t) 

(36) 

(37) 

(38) 

(39) 


Substituting  (25)  into  (19)  and  using  the  above  identities, 
one  obtains  the  equations  of  motion  for  the  coefficients  of 
the  perturbation  expansion: 

Tr  ^  Aw  =  [ - 2gl  +  2yA *]Aw  +  TR SW(T)  (40) 


Fig.  I.  Timing  fluctuations  as  computed  from  our  theory  and  by  a  (nu¬ 
merical)  integration  that  includes  back  coupling  of  continuum;  the  timing 
is  normalized  to  the  noise  source.  Solid  line  the  analytic  result,  dashed  line 
computer  result;  g ,  =  0. 

From  the  preceding  equations  one  draws  the  following 
conclusions:  The  energy  will  grow  exponentially,  unless 
the  gain  saturation  is  sufficient  to  suppress  the  APM  ac¬ 
tion,  as  indicated  by  (40).  An  energy  change  couples  to 
the  phase  evolution,  because  the  change  affects  the  Ken- 
phase  shift,  according  to  (41).  Frequency  deviations  damp 
out,  because  the  gain  of  a  pulse  whose  center  frequency 
is  off  the  gain  line  center  is  nonuniform  across  the  spec¬ 
trum  and  pushes  the  spectrum  back  to  line  center.  Finally, 
a  change  of  frequency  and  a  change  of  energy  both  cause 
a  change  of  group  velocity,  the  first  because  of  GVD,  the 
second  since  by  the  Kramers-Kronig  relation  the  gain 
profile  is  associated  with  an  index  profile  that  changes 
along  with  the  gain. 

As  mentioned  earlier,  the  coupling  back  from  the  con¬ 
tinuum  via  gain  dispersion  and  APM  or  KLM  has  been 
neglected  causing  an  error  of  order  fi1.  Fig.  1  shows  a 
comparison  of  the  prediction  of  our  theory  with  the  exact 
numerical  solution  of  the  perturbation  equation  (16)  for 
assumed  independent  noise  sources  Sj.  The  numerical 
method  of  solution  has  been  discussed  elsewhere  [15].  The 
ordinate  is,  of  course,  proportional  to  the  intensity  of  the 
noise  source  and  is  normalized  to  it.  It  should  be  noted 
that  the  agreement  is  indeed  excellent. 


V.  The  Noise  Sources 

We  consider  separately  the  noise  introduced  by  the 
fluctuations  of  particular  physical  quantities. 


Tr  ~  A6  =  -6 A20  —  +  TRSe(T)  (41) 

dT  w„ 

=  -\^*p  +  Wr>  (42) 

Tr  ^  =  -2| D  |  Ap  -  -f-  —  +  TrS,(T).  (43) 
dT  n,  w„ 

The  noise  sources  are  defined  by 

Sj(T)  =  \  {  dt  [f*(t)S(t,  T)  +  T)\  (44) 


Gain  Fluctuations 

As  one  can  see  from  the  master  equation  (1).  a  (long 
term)  fluctuation  A g  ( T )  of  the  gain  introduces  an  equiv¬ 
alent  noise  source: 


SaA  T) 


*g(T) 

Tr 


1  d\ 

-Qjt)^' 


T). 


(45) 


When  the  laser  frequency  depends  on  the  gain,  which  is 
the  case,  for  example,  when  different  spectral  shapes  of 
the  gain  and  the  loss  force  the  laser  to  operate  away  from 
the  maximum  of  the  gain  curve  gain  fluctuations  give 
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also  the  frequency  pulling  term  (see  Fig.  2): 


•Wri  T)  -  j 


Ag(T) 


47*  («o  ~  dt 


-as«,T).  (46) 


This  term  does  not  follow  directly  from  (1),  but  is  ob¬ 
tained  by  extension  of  the  gain  operator  to  carrier  fre¬ 
quencies  u0  off  line  center. 


Length  and  Refractive  Index  Fluctuations 

Fluctuations  of  the  cavity  length  give  rise  to  two  noise 
contributions: 


SAL(t,  T) 


*L(T) 

Tr 


1  d 
n  +  —  — 
vgdt 


as{t,  T ) 


(47) 


where  n  is  the  effective  index  assigned  to  the  total  length 
L. 

The  first  term  in  parentheses  is  the  frequency  pulling 
effect  of  a  cavity  of  varying  length.  The  second  term  gives 
the  timing  jitter  caused  by  a  cavity  of  varying  length. 
Cavity  geometry  fluctuations  in  a  KLM  system  may  also 
cause  fluctuations  of  the  effective  saturable  absorber  pa¬ 
rameter  7,  if  the  cavity  is  so  laid  out  as  to  enhance  it  by 
operating  near  the  stability  boundary  of  the  cavity  mode. 
If  L  is  again  a  typical  cavity  dimension,  A Ldy /dL  is  an 
effective  noise  term  in  the  master  equation. 

Fluctuations  of  the  refractive  index  produce  noise  terms 
that  are  analogous  to  those  due  to  cavity  length  fluctua¬ 
tions,  with  A L  replaced  by  L(An/n). 


Spontaneous  Emission  Fluctuations 

Beside  the  classical  noise  sources  considered  above, 
spontaneous  emission  fluctuations  are  always  present  into 
the  cavity  for  reason  of  quantum  mechanical  consistency. 
In  a  semiclassical  approach,  the  noise  introduced  by  spon¬ 
taneous  emission  in  a  laser  is  taken  into  account  introduc¬ 
ing  a  white  noise  source  with  correlation  [16],  [17] 

<S,„( 7-,  t)S*na\  (')>  =  0^hv8(T  -  T')8 (t  -  f) 


where  0  is  the  enhancement  factor  due  to  the  incomplete 
inversion  of  the  medium. 


Noise  Sources  in  the  Equations  of  Motion 
By  projecting  via  the  adjoint  set  of  functions,  we  get 
Ag(T) 


SJT)  =  2 


wa  +  Swg„(T) 


(49) 


5,(7-)  = 


SAT)  = 


vJr 


AL(T)  +  L 


A/r(D 


n 


Ag(T) 


3rR(w0  -  Mw|r 


+  Se,,AT)  (50) 
2  +  S^AT)  (51) 


Fig.  2.  Gain  and  loss  profiles  with  different  center  frequencies  for  which 
gain  fluctuations  cause  frequency  fluctuations. 


S,(T) 


*g(T) 


1 

7>* 


AL(T)  +  L 


An(T) 

n 


+  S,,qn(T)  (52) 

where  the  terms  Si  qn(T),  due  to  the  quantum  noise,  are 
independent  white  noise  terms 


< S,,„(T)SMn(r)>  =  D,  qn8j  hb(T  -  T')  (53) 

with  diffusion  constants 


2  g 


D*.qn  =  4wo0  —  hv 
*  n 


4  ( .  rL\  .  2g 


D-  =  ^{'+T2)efA 


(54) 

(55) 


D, 


p.qn 


D,. 


qn 


2  4 

3 w„T2 


At2 


3  w„ 


(56) 

(57) 


For  white  noise  sources  it  is  convenient  to  write  the  cor¬ 
relation  functions  as 


(SATySjlT’))  =  Dj  j5(T  —  T').  (58) 

The  diffusion  constants  Dt  J  form  a  positive  semidefinite 
matrix.  The  assumption  of  white  noise  simplifies  the  an¬ 
alytic  treatment.  Anyway,  one  should  be  aware  of  the  fact 
that,  if  quantum  noise  is  white,  the  classical  noise  sources 
in  general  are  not.  The  validity  of  the  results  found  as¬ 
suming  classical  white  noise  is  then  restricted  to  the  fre¬ 
quency  region  inside  the  actual  bandwidth  of  the  noise 
sources.  In  the  ensuing  analysis,  the  diffusion  coefficients 
will  be  assumed  as  given,  and  not  necessarily  as  those  of 
quantum  nature  (54),  (57). 


VI.  The  Noise  Spectra  and  the  Correlation 
Functions 

The  equations  of  motion  (40)-(43)  are  of  two  types. 
One  type  is  a  relaxation  equation  driven  by  a  noise  source, 
the  other  has  no  relaxation  time,  and  is  driven  by  a  noise 
source  as  well  as  one  of  the  perturbation  variables.  In  the 
former  case,  if  the  relaxation  times  are  relatively  long 
compared  with  the  inverse  bandwidth  of  the  noise  sources. 
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the  spectrum  of  the  noise  sources  can  be  treated  as  white 
and  the  resulting  fluctuations  have  a  Lorentzian  spectrum. 
Let  us  define  the  Fourier  transform  pair  as 


f(T)  =  >!¥„  J  ^  e^Q)  (59) 


where  the  infinitely  long  normalization  time  Tc  has  been 
introduced  to  avoid  the  divergence  of  the  fluctuation  spec¬ 
tra  which  is  present  with  the  usual  definition  of  Fourier 
transform.  This  divergence  arises  from  the  fact  that  delta 
correlation  in  time  domain  corresponds  to  delta  correla¬ 
tion  in  frequency  domain.  The  inverse  of  the  time  T„  gain 
the  “width”  of  the  delta  functions  in  the  frequency  do¬ 
main,  and  renormalizes  their  amplitude. 

Let  us  define  the  two  relaxation  times: 

-  =  (2g,  -  2yAl)^r  (60) 

Tw  Tr 


1  =  i  8 

r„  3  TRSiy 


The  spectra  of  the  energy  fluctuations  are  then 


<|6w(n)|J) 


<|S„(«)|2> 


(61) 


(62) 


and  of  the  frequency  fluctuations 


<|Ap(ft)|2> 


<|S„(ft)|2> 


(63) 


These  spectra  couple  to  the  phase  and  timing  fluctuations. 
If  the  noise  source  spectra  are  treated  as  white,  these  latter 
spectra  are  of  the  shape: 


which  can  be  expanded  into  their  respective  singularities. 
The  leading  singularity,  which  gives  the  asymptotic  be¬ 
havior  for  large  T ,  has  a  spectrum  of  the  form  1  /fl2,  cor¬ 
responds  to  a  random  walk.  Concentrating  on  this  domi¬ 
nant  term,  one  finds  that  the  phase  and  timing  mean  square 
fluctuations  are  of  the  form: 

<| At  (7  +  T„ )  -  At(Ta)\2)  =  D,T  (65) 


with 


D,  =  <|Sr(n  =  0)  | 2  > 


+ 


4Z)2(l5p(n  =  0)| 2 ) 

Tr/ Tp 


and 


<|A0(7  +  T„)  -  A6(Tn) |2>  =  DeT 


(66) 

(67) 


with 


P,  -  -  0)1’)  +  <8/l“)!<l,S*,(^  7  °>|;>  (68) 


where  the  coefficients  multiplying  T  play  the  role  of  dif¬ 
fusion  coefficients. 

The  complete  expressions  for  the  correlation  functions 
of  timing  and  phase  fluctuations  is  more  involved.  They 
are  anyway  easily  calculated  from  the  fluctuation  spectra. 

The  complete  expression  for  the  timing  jitter  spectrum 
is 


<|A,(fl)|2>=^^«>  L  mm2) 


ft2  (  ft2  +  ~ 


n 


and  for  the  phase  noise  spectrum 

(S/42)2  <|S„.(ft)|2> 


<|A0(n)|2>  = 


n2  (n2  +  V 

V  rij 


n2 

(69) 

<|S0(O)l2) 

Q2 


(70) 


The  timing  mean  square  fluctuations  can  be  easily  ob¬ 
tained  by  using  the  relationship 

<|At(7  +  T„)  -  Af(7„)|2> 

=  ~-[  < | Af (12) (e /or  -  1)|2>  <Kl 

lit  J  -Oo 

=  -  (  <  I Ar  (12)) 2  >  11  -  cos  ((IT)]  d(l.  (71) 

ir  Jo 

Assuming  white  noise  spectra  for  the  frequency  and  tim¬ 
ing  noise  sources,  we  get 

<|At(T  +  T„)  -  At(Ta)\2) 

=  ~  Dp,/p  (-  -  1  +  e~T/A  +  D,  ,T.  (72) 

lR  Jp  / 

It  is  interesting  to  calculate  the  two  asymptotic  limits  of 
(72).  Neglecting  the  contribution  of  D,  „  for  T/rp  «  1 
we  have 

.  4 D2  Dp.pTPT2 

<1  A/(7  +  T„)  -  At(T„)Y >  *  -r  ■■  (73) 

*  R  *■ 

while  for  T/tp  »  1 

,  4  D2  , 

<|At(7  +  T„)  -  At(7„)|2>  *  -=r  Dp  pt2T  (74) 

1  R 

The  variance  of  the  timing  fluctuations  is  unbounded,  as 
expected  from  the  absence  of  any  time  reference  in  pas¬ 
sively  modelocked  lasers.  It  grows  for  small  times  qua- 
dratically  and  for  large  times  linearly  with  T.  The  quadratic 
growth  with  time  arises  from  the  i  /  I)4  dependence  for 
n  »  i/v  while  the  linear  dependence  corresponds  to 
the  1  /  ft2  spectrum  for  (l  «  1  /rp. 
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It  is  worth  noting  that  for  negligible  initial  frequency 
fluctuations  a  direct  integration  in  the  time  domain  of  (42) 
and  (43)  gives,  in  the  asymptotic  limit  T  «  tp 

AD2 

<|A t(T  +  T0)  -  Ar(70)|2>  *  -2-  DppT\  (75) 

1  R 

different  from  the  previous  result.  This  T 3  dependence  is 
closely  related  to  the  cubic  dependence  of  timing  fluctua¬ 
tions  on  propagation  distance  in  long  haul  soliton  com¬ 
munication  systems  known  as  the  Gordon-Haus  effect 
[16].  The  different  asymptotic  behavior  of  (73)  is  due  to 
the  different  natures  of  the  two  processes:  in  a  laser  we 
are  dealing  with  steady-state  properties,  that  do  not  exist 
in  the  propagation  of  a  soliton  in  a  fiber,  by  its  nature  a 
transient  process. 

The  phase  mean  square  fluctuations  can  be  evaluated  in 
a  similar  way,  obtaining 

<|A0(F  +  Ta)  -  A8(T0) |2> 


(25  A])2 


—  - 1 
T*. 


(76) 

Other  spectra  and  correlation  functions  are  of  interest.  We 
will  assume  for  the  sake  of  simplicity  that  the  noise 
sources  are  white.  From  the  spectrum  of  energy  fluctua¬ 
tions,  we  get  the  correlation  function  of  the  energy  fluc¬ 
tuations 

<w(7  +  T0)w(T0))  =  ~  J  „  **r<M0)|2>  M 


=  yA*.  H,e',r|/I“.  (77) 

The  cross-correlation  spectrum  of  time  and  frequency 
fluctuations  reads 


<A/*(«)Aw(fl)>  = 


_ -2\D\DP,„ 

-jU(-jQ  +  +  1/rJ 


and  the  cross-correlation  function 
<[A/(7  +  Ta )  -  Ar(70)]Aw(r0)> 


.±r 

2t  J-o 


_ 


l)<Ar*(fl)Aw(Q)>  dSl 


=  -2  sign  (T)\D\D.„ 


Tp  “l"  Tw 


t ( 1  -  e-|r|/T) 


where  t  is  equal  to  t_  for  T  >  0  and  to  tw  for  T  <  0. 


VII.  The  Measurement  of  Noise 
The  most  direct  way  of  measuring  the  different  noise 
spectra  would  be  to  construct  a  homodyne  detector  with 
a  pulsed  local  oscillator  having  one  of  the  adjoint  func¬ 


tions  /  as  the  pulse  shape.  This  pulse  would  be  derived 
from  the  mode-locked  pulse  and  delayed  by  the  time  T. 
Successive  samples  of  the  square  of  the  detected  output 
would  be  collected.  Thus  one  would  arrive  at  the  mean 
square  fluctuations  of  the  four  pulse  parameters.  In  some 
laser  systems,  the  delays  would  have  to  be  considerable, 
since  the  relaxation  times  of  some  solid-state  lasers  are 
many  microseconds  long.  Delay  fiber  “lines”  would  dis¬ 
tort  the  pulses  and  the  experiment  would  not  be  all  that 
simple.  An  alternative  way  is  based  on  the  detection  of 
the  low  frequency  spectrum  of  the  detected  pulse  [21]. 
Let  us  concentrate  on  this  approach. 

The  periodic  pulse  train  has  the  time  dependence: 

a(r,  T)  =  E  A0  sech  At  ~  mT^j 

(80) 

The  field  amplitude  a(t,  T)  is  a  function  of  two  time  vari¬ 
ables,  t  and  T\  the  former  is  the  fast  time  variable  in  terms 
of  which  the  pulse  envelope  is  expressed;  the  iatter  is  the 
time  variable  of  the  slow  time  scale  in  terms  of  which  the 
noise  is  expressed.  The  dependence  on  T  is  implicit  in  the 
time  dependences  of  the  perturbations  of  energy,  Aw(7), 
of  frequency,  A p(T),  of  phase,  A 9(T),  and  of  timing, 
A t(T).  The  pulse  repeats  every  round-trip  time  TR.  Since 
the  laser  relaxation  rates  are  slow  compared  with  one 
round-trip  time,  the  fluctuations  produce  narrow-band 
spectra  around  each  of  the  Fourier  components  of  the  pe¬ 
riodic  pulse  train  separated  by  0o  =  2-k/Tr.  One  may 
decompose  the  function  a(t,  T)  into  a  sum  of  Fourier 
components: 

a(t,  T)  =  x  —  E  sech  ^  tvSX0 
•r  «  L 2  J 

.  eJ(nO<,  +  a0- &pV  g-jnOoM 

When  the  pulse  train  impinges  on  a  detector,  the  kth 
Fourier  component  of  the  detector  current  is  produced  by 
the  beat  of  the  nth  Fourier  component  of  the  amplitude 
with  the  (n  +  k)th  Fourier  component: 

l,('-  n,‘  *  ”  i  (’  if) 

•  sech  y  [(m  +  k)Qa]  e’™*  e',ka (82) 

where  i}  is  the  quantum  efficiency  of  the  detector  and  e  is 
the  electron  charge.  We  note  that  the  phase  fluctuations 
A  6  and  the  frequency  shift  Ap  do  not  appear  in  the  Four¬ 
ier  components  of  the  current.  In  order  to  measure  them, 
one  needs  a  different  experiment.  We  shall  return  to  the 
discussion  of  such  a  measurement  further  on. 

Since  the  pulsetrain  contains  very  many  Fourier  com¬ 
ponents  (1000  or  more),  one  may  expand  the  expressions 
to  first  order  in  krft0.  Further,  one  notes  that  the  pulse- 
width  itself  is  subject  to  fluctuations,  along  with  the  am¬ 
plitude  A„.  If  the  pulse  is  soliton  like,  its  area  is  invariant 
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and 

A  (wt)  =  A(2/12t2)  =  0 

so  chat 

A  w  At 

w  T 


(83) 

(84) 


The  amplitude  of  the  kth  Fourier  component  of  the  current 
becomes: 

•  1  -  ~  Tkti0  tanh  rmSl^ 

■  1  -  2  y  ^  Tm®°  j  tanh 

-  y  (|  tanh  (85) 

One  may  replace  the  sum  over  m  by  an  integral .  Only 
the  symmetric  functions  of  m  contribute.  The  result  is 

ik(T)  =  91  1+/*—  e~J^  (86) 

w 


where 


Consistently  with  the  first-order  expansion  in  xkrf]0/2  we 
will  neglect  the  second  term  in  the  definition  of  fk  by  set¬ 
ting  /*  =  1 .  If  the  energy  and  timing  fluctuations  are  in¬ 
dependent,  the  correlation  function  of  the  current  is  the 
sum  of  the  correlation  functions  of  exp  ( —j  kft„A  t )  and  of 
At,  and  so  is  the  spectrum.  If  the  timing  fluctuations  go 
through  a  random  walk,  the  spectrum  of  exp  (—j  kfl„Ar) 
is  Lorentzian.  The  spectrum  of  At  is  added  to  it.  The 
power  of  the  Fourier  components  of  the  former  increases 
with  k2,  the  latter  has  a  constant  contribution  plus  one  that 
increases  with  k2  and  another  that  increases  with  k4.  Thus 
the  two  fluctuations  A w(7)  (or  A t(7))  and  At (7)  con¬ 
tribute  to  the  current  spectrum.  They  are  independently 
measurable.  At (7)  experiences  a  random  walk  with  the 
diffusion  constant  D, .  Hence  it  produces  a  linewidth  con¬ 
tribution  to  the  kth  Fourier  component  that  is 

Afl  =  ?D,/2  (89) 

The  latter  has  a  contribution  that  is  independent  of  k. 
Thus  the  two  fluctuations  Aw (7)  (or  At(7))  and  Ar(7) 
contribute  to  the  current  spectrum.  They  are  indepen¬ 
dently  measurable  on  the  basis  of  their  different  k  depen¬ 
dence. 

When  A  /  and  A  t  are  correlated,  and/or  the  timing  fluc¬ 
tuations  are  not  those  of  a  simple  random  walk,  the  spec¬ 


trum  is  more  complicated.  Let  us  define  the  spectrum  of 
the  amplitude  of  kth  Fourier  component  of  the  current 
fluctuations  as 


=  (  e~JQTCk(T)  dT 

—  03 


where 


C*(7)  =  <4(7  +  T0)it(T0))  .  (91) 

Since  C*(7)  =  C*(—  7),  the  spectrum  can  also  be  writ¬ 
ten  as 


/*(«)  =  2  Re 


f  e"'QrC* 

Jo 


(7)  dT.  (92) 


The  correlation  function  follows  from  (86) 


C*(7)  =  9l2  (  1  +  -  [A  w(7  +  T0)  +  A  w*(70)] 


+  —  [Aw(7+  70)Aw*(70)] 
w 


.  e-jkQA&HT 


+  t„)-a/(7„)]^ 


For  Gaussian  processes  the  average  can  be  evaluated  [20] 
and  the  result  is 

Ck(T)  =  3l2^l  +  y  <A w(T  +  T0)Aw(T„)> 


- f  <[A t(T  +  T„)  -  Ar(70)]A w(7)> 

w 

•  <[Ar(7  +  T0)  -  A/(70)]Aw(70)> 


-j  —  [<[A/(7  +  Ta)  -  At(Ta)]A  w(7)> 
w 


+  <[A/(7+  Ta)  -  A/(70)]Aw(70)>] 


exp  --k2fl2<[A/(7+  7„)  -Ar(7„)]2 >.  (94) 


Substituting  in  this  expression  the  values  of  the  correla¬ 
tion  functions,  and  taking  into  account  that 

<[Ar(7  +  T0)  -  A/(70)]Aw(7)> 

=  — <[A/(  — 7+  T0)  -  At(T0)]Aw(T0))  (95) 

we  get,  for  7  >  0  and  for  white  noise  sources 

Ck(T)  =  91 2  {1  +  fl,e“(r/r*)  -  a\&Tprwll  -  e',r/r'’ ) 

•  (1  -  e~{~T,Tr)  +  ka2[TH.(l  -  e_<r/T*’) 

+  t„(1  -  e-'W))} 

•  exp  -k2a3  (~  ~  1  +  e~a/Tp)^  (96) 
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where 


1  Tw  „ 

C|  '  ^  2  ^w  w 

c2  =  2\D  |  ft0  £)„ 

W/R  7>  +  Tw 

2D 2 

ci  —  —2  Dp,p®oTp- 

‘  R 


Dw,„Dpp  >  Dwp  (100) 

for  the  positive  semidefiniteness  of  the  correlation  matrix, 
the  coefficients  c,  have  to  meet  the  condition 

->  -Z- (Tp  +  Tw?c\  (101) 

"  *  H' 

where  the  equality  sign  holds  for  perfectly  correlated  noise 
sources,  when  energy  and  frequency  fluctuations  are 
driven  by  the  same  noise  sources  (e.g.,  gain  fluctuations). 

The  spectrum  can  be  easily  evaluated  analytically,  since 
the  only  integral  which  is  involved  can  be  easily  per¬ 
formed  with  the  substitution  e  x  -*  y, 

(  exp  (-ax  -  be~x)  dx  =  b~“y(a,  b)  (102) 


Frequency  oflttt  (Hz) 

(a) 


V 

Frequency  offMt  (Hz) 

(b) 


where 


7  (a,  b)  =  f  e  'ta  1  dt 
Jo 


is  the  incomplete  gamma  function.  If  we  define 

<Ma,  ft)  =  Tpek2c^ai)(jari'-a-k2ci) 

■  7[-(jflTp  -  a  -  ^c3),  A^c3]  (104) 

the  spectrum  of  the  Irth  Fourier  component  becomes 


/*(0)  =  291  Re  j^(0,  0)  +  c, <t>  (y,  0^ 

-  c22kl7pTw  <t>( o,  0)  -  <t>  nj  -  <t>( i, 

+  *(^+lt0). 


+  jkc2  (tp  +  7v)<f>(0,  0)  -  7.0(1,  0) 


This  expression  has  been  used  to  plot  the  curves  shown 
in  Fig.  3.  When  intensity  and  timing  fluctuations  are  un¬ 
correlated,  for  c3  =  0,  the  spectrum  is  symmetric.  Asym¬ 
metry  appears  when  c3  *  0  when  intensity  and  phase  fluc¬ 
tuations  are  correlated,  as  shown  in  Fig.  4. 


-510*  0  510* 

Frequency  oftet  (Hz) 

(c) 

Fig.  3.  The  spectra  of  different  current  Fourier  components  for  uncorre- 
lated  intensity  and  timing  fluctuations.  The  values  of  the  parameters  are: 
tw  =  0.5  ms.  rf  -  8.4  ms,  c,  -  25  x  I0~4  (5%  intensity  fluctuations), 
c,  =  0,  c,  =  0.0073. 

It  is  interesting  to  give  a  simplified  expression  of  i*( fl) 
for  larre  0.  Large  0  corresponds  to  small  T,  for  which 
the  exponential  in  (86)  can  be  expanded  to  first  order  ob¬ 
taining 

ik(T)  *  91  1  +  —  Aw(D  1  -  jkQ0At(T)  . 

I  w  J 

(106) 

From  this  expression  the  spectrum  is  directly  evaluated 

<|i4(0)|J>  *  91)2x6(0)  +  —2  <|Aw(0)|2> 

+  *r2n2  <  j  Ar(0)|2  > 

+  jkno[<  A/*(fl)Aw(0)> 


-  <  Ar(Q)Aw*(Q)>] 


Von  der  Linde  was  the  first  to  propose  the  use  of  (107) 
to  measure  the  energy  fluctuations  and  the  timing  jitter 
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Frequency  oflMt  (Hz) 

(a) 


Frequency  otliet  (Hz) 

fc) 


Fig  4  The  spectra  of  differ  it  ...  Trent  Fourier  components  for  correlated 
intensity  and  timing  fluctuations  The  values  of  the  parameters  are  the  same 
as  in  Fig.  3.  and  r*  =  V4< ,  c,r. /lr„(rr  +  r.)’|  |see  (10D). 

spectra  in  mode-locked  lasers.  He  considered  the  case  of 
uncorrelated  energy  and  timing  fluctuations,  when  the 
quantity  within  square  brackets  in  (107)  is  zero.  He  noted 
that  the  k  th  Fourier  component  of  the  detected  current  is 
a  Dirac  delta  function,  corresponding  to  the  noiseless 
laser,  plus  the  sum  of  the  energy  fluctuation  spectrum  and 
of  k2  times  the  timing  jitter  spectrum.  The  contributions 
of  energy  fluctuations  and  timing  jitter  can  be  resolved  by 
using  their  different  k  dependences.  In  the  original  paper, 
the  theory  was  developed  mainly  for  active  modelocking 
where,  due  to  the  active  modulation,  timing  fluctuations 
are  a  stationary  process.  Later  on,  this  theory  has  also 
been  applied  to  the  study  of  the  time  jitter  of  passive 
mode-locked  systems.  Since  in  passive  mode-locking  the 
timing  fluctuations  are  not  a  stationary  process,  the  von 
der  Linde's  method  has  to  be  used  carefully.  To  see  the 
differences,  it  is  useful  to  consider  the  case  in  which  the 
timing  jitter  is  a  stationary  process,  by  introducing  a  re¬ 
storing  force  in  the  equation  of  motion  for  &t(T).  This 
accounts  for  the  introduction  of  a  modulator,  locked  to  an 


external  oscillator,  into  the  laser  cavity  as  proposed  in 
[22]  to  achieve  self  starting  and  reduction  of  the  timing 
jitter.  Equation  (43)  is  modified  by  adding  the  term 
-A t/r,  to  the  RHS,  where  r,  is  the  decay  time  of  the 
timing  fluctuations.  The  spectrum  of  the  timing  jitter  be- 


<|  Ar(fl)|2>  =  -r 


<IVft)l2> 


0*  +  "2 


)("*) 


n2  +  \ 

The  1  /fl2  divergence  typical  of  random  walk  present  in 
(69)  is  removed  as  expected  because  timing  jitter  is  now 
a  stationary  pr.  cess.  The  timing  jitter  correlation  function 
becomes 

<|A t(T+  T0)  -  Ar(7'„)|2 > 


t-2  Dp.p  2 
I  R  r  t 


r,  ~  rp 


[7,(1  _  ,-iri/r() 


-  7.(1  -  e-|r|/")]  +  D,,t,(1 


_  e~\T\/r,) 


The  expression  for  the  correlation  function  of  the  kth 
Fourier  component  of  the  photocurrent,  in  the  case  of  un- 
correlated  time  and  energy  fluctuations,  is 


Ct(T)  =  91 2  p  ^  <  A wfT  4  T„)Aw(T0))  | 

•  exp  -^k2n2<[Ar(T+  T„)  -  A r(T0)]2  >. 


The  correlation  function  of  the  energy  fluctuations  is  given 
by  (77).  The  spectrum  of  the  kth  Fourier  component  is 
given  by  the  Fourier  transform  of  Ct(T).  It  is  interesting 
to  see  what  happens  for  |T|  -*  oo:  the  timing  fluctuations 
given  by  (109)  approach  asymptotically  a  constant  value. 
This  means  that  the  fluctuation  spectrum  at  0  =  0  is  di¬ 
vergent.  Let  us  examine  the  nature  of  this  divergence  by 
considering,  for  simplicity,  the  case  of  r,  »  tp.  In  this 
limit,  the  fcth  component  of  the  photocurrent  spectrum  be- 


<!«*( 0)1 2  > 


(1  + 


-<iri/T,), 


•  exp  [-^^7,(1  -  f-,|rl/T'>))  dT  (111) 

where  dt  and  d}  are  positive  constants  proportional  to  the 
noise  diffusion  constants: 

,  r„  Dw 

,'1 '  '2?r 


2Q;,D1Dp,pt; 
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By  expanding  exp  [d3k2T,e  (|r|/T,)]  inio  a  series  and  inte¬ 
grating  we  get 


r 


<|i*(0)|  >  =  * 


2  \  —  e-d*1'! 


2x6(12)  +  — - — 

r“  n2  1 


Qz  +  — 
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+  s 


(d3  fcV 


n  =  1  n] 


2  2  n 

—  +  — 
Tu-  r, 


Bz  +  I  —  +  — 


fi  4*  - n 


(112) 


For  12  =  0  the  series  is  convergent,  and  the  only  term  left 
is  the  Dirac  delta  function.  We  have  found  that  the  diver¬ 
gence  of  the  current  spectrum  at  12  =  0  is  due  to  a  delta 
function.  Equation  (107)  is  obtained  by  keeping  in  the 
exact  expression  only  the  terms  of  first  order  in  d,  and  d3. 
Summarizing,  the  *th  order  of  the  low  frequency  noise 
spectrum  of  the  photocurrent  is  a  delta  function  at  12  =  0 
plus  a  broad-band  noise  spectrum  which,  at  first-order,  is 
the  sum  of  the  energy  and  k2  times  the  timing  jitter  spec¬ 
tra.  The  first  order  expansion  of  the  exponential  that  leads 
to  (107)  is  then  justified,  because  the  terms  left  over  in 
the  expansion  are  small  for  all  values  of  12.  In  Fig.  5  the 
result  of  the  approximate  and  exact  expressions  are  com¬ 
pared.  The  solid  line  is  the  plot  of 


i*  (12)  =  e'^'  2t,  Re 


(-d3^T,)'jaT'7(y'12r„  -  d3k?7,) 


.  y  (ja,,  +  1.  _ 


(113) 


which  is  equivalent  to  the  sum  of  all  the  terms  of  (112) 
and  has  indeed  a  delta  function  singularity  in  the  origin, 
while  the  dashed  line  is  the  plot  of  the  only  first-order 
terms  in  and  d3.  The  approximate  and  exact  analyses 
give  indistinguishable  results. 

The  general  structure  of  the  photocurrent  spectra  with 
a  Dirac  delta  function  and  a  broad-band  background  is 
common  to  all  the  cases  in  which  the  presence  of  a  re¬ 
storing  force  in  the  timing  fluctuations  makes  the  timing 
jitter  a  stationary  process  like,  e.g.,  the  active  mode  lock¬ 
ing. 

In  the  case  of  pure  passive  mode-locking  we  have  a 
different  picture.  In  this  case,  the  linear  divergence  of  the 
timing  fluctuations  for  long  |  7j  due  to  the  absence  of  any 
restoring  force  prevents  the  Fourier  transform  of  Ck(T)  to 
diverge  for  12  =  0.  This  means  that  now  the  Dirac  delta 


Fig.  5.  The  spectra  of  different  current  Fourier  components  in  presence  of 
an  active  mode-locking  restoring  force.  The  values  of  the  parameters  are: 
r„  =  0.5  ms.  r,  =  1  ms,  if,  =  25  x  1(T4,  d ,  =  0.0073. 


function  at  12  =  0  of  the  approximate  expression  (107)  is 
nonphysical,  and  it  is  present  only  because  we  have  per¬ 
formed  a  first-order  expansion  of  the  exponential  for  times 
where  this  approximation  is  not  allowed.  Equation  (107) 
is  a  good  approximation  of  the  actual  photocurrent  spec¬ 
trum  only  down  to  a  cutoff  frequency  where  the  variance 
of  the  timing  fluctuations  times  it2!!*  becomes  of  the  order 
of  1.  This  means  that,  in  principle,  the  timing  jitter  spec¬ 
trum  at  very  low  frequency  cannot  be  directly  measured 
by  the  method  proposed  by  von  der  Linde  which  is  based 
on  the  validity  of  the  linearized  equation  (107).  The  ap¬ 
proximate  and  the  exact  expressions  of  the  photocurrent 
spectrum,  calculated  assuming  zero  energy  fluctuations, 
are  compared  in  Figs.  6  and  7  assuming  different  values 
for  the  diffusion  constant  for  the  frequency.  In  Fig.  6  the 
values  of  the  parameters  have  been  chosen  to  fit  the  results 
of  [22].  One  sees  that  the  spectrum  rises  proportionally  to 
k 2,  and  that  the  comparison  between  approximate  and  ex¬ 
act  expressions  is  indeed  excellent  down  to  100  Hz.  Since 
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Fig.  6.  The  spectrum  of  the  current  fluctuations  of  the  first  and  20th  Fou¬ 
rier  component  of  the  current.  The  values  of  the  parameters,  chosen  to  fit 
the  results  of  (22J,  are  tp  =  8.4  ms,  c,  =  0.  Cj  =  0  (zero  intensity  fluc¬ 
tuations),  r,  =  0.0073.  Dashed  approximate,  solid  exact. 


Fig.  8.  The  spectrum  of  the  current  fluctuations  for  the  first  and  20th  Fou¬ 
rier  component  for  the  same  values  of  Fig.  7  and  t_  =  0.5  ms,  c2  =  25  x 
!0'4  (5%  intensity  fluctuations).  Dashed  approximate,  solid  exact. 


Fig.  7.  The  spectrum  of  the  current  fluctuations  of  the  first  and  20th  Fou¬ 
rier  component  for  higher  values  of  the  frequency  diffusion  constant.  The 
values  of  the  parameters  are:  rp  =  8.4  ms.  c,  =  0.  c,  =  0.073. 


the  cutoff  frequency  below  which  the  approximate  expres¬ 
sion  fails  is  proportional  to  k2,  we  see  a  large  deviation 
from  the  exact  theory  only  when  we  consider  very  high 
harmonics. 

When  timing  jitter  increases,  the  approximate  and  exact 
analyses  give  different  results  for  low  harmonic  numbers. 
This  is  shown  in  Fig.  7,  where  the  plots  of  Fig.  6  are 
repeated  for  a  ten  times  larger  value  of  the  frequency  dif¬ 
fusion  constant.  In  this  last  case  the  exact  and  the  approx 
imate  analyses  give  different  results  even  for  n  =  20. 
Consequently,  the  shape  of  the  timing  jitter  spectrum 
(given  by  the  approximate  analysis)  cannot  be  obtained 
by  the  analysis  of  the  high-order  components  of  the  pho¬ 
tocurrent  spectrum.  On  the  other  hand,  the  lower  order 
components  are  affected  by  the  intensity  noise  as  shown 
in  Fig.  8  where  the  case  of  Fig.  7  has  been  plotted  by 
adding  an  intensity  noise  of  5  % . 

Phase  and  frequency  fluctuations  need  a  different  way 
of  detection.  A  cascade  of  Fabry-Perot  interferometers 
can  isolate  the  fundamental  n  =  0  of  the  pulse  spectrum 

ajt,  T)  =  *  —  (114) 

When  the  diffusion  of  the  phase  can  be  approximated  to  a 
random  walk,  it  causes  a  Lorentzian  shape  of  the  fre¬ 
quency  spectrum  whose  width  is  equal  to  the  diffusion 
constant  of  the  phase  [23].  When  this  approximation  can¬ 


not  be  performed,  its  complete  expression  is 
Sm((1)  =  I  e-JQTe~<i/2)<WT  +  ra)-&t(Tani)  dj 

J  — OD 

=  2  Re  ~  *.),  6.U 

(115) 


where 

*>  (H6) 
*woJ  K 

The  spectrum  of  is  Gaussian  and  its  expression  is 

(117) 

The  net  spectrum  is  the  convolution  of  the  two  spectra. 


VIII.  Experimental  Results 
Let  us  examine  briefly  the  existing  experimental  data 
on  timing  jitter  of  passive  mode-locked  lasers.  Many  ex¬ 
perimental  results  have  been  published  so  far  on  timing 
jitter  of  colliding  pulse  mode-locked  (CPM)  lasers  [24], 
[25]  and  of  color-center  lasers  [26].  Even  though  our 
model  applies  to  APM  and  KLM  lasers,  all  reported  ex¬ 
perimental  configurations  involve  intracavity  negative 
dispersion  to  compensate  for  the  self-phase  modulation. 
The  mechanism  that  produces  the  coupling  between  fre¬ 
quency  fluctuations  and  timing  fluctuations  (the  depen¬ 
dence  on  frequency  of  the  group  velocity)  is  expected  to 
be  active  also  in  these  cases.  Generally,  it  is  difficult  to 
obtain,  from  the  published  experiments,  all  the  informa¬ 
tion  required  on  the  numerical  values  for  the  parameters 
of  the  theory  to  make  quantitative  comparisons.  How¬ 
ever,  the  very  specific  qualitative  predictions  of  the  theory 
can  be  compared  with  available  experimental  informa¬ 
tion.  In  ail  the  reported  timing  jitter  spectra,  when  no  par¬ 
ticular  arrangement  is  used  to  reduce  the  timing  jitter,  the 
power  spectra  drop  40  dB /decade,  which  corresponds  to 
the  expected  1  /(T  behavior.  More  interesting  is  the  jitter 
measurement  in  [22],  where  the  investigated  laser  closely 
corresponds  to  our  model.  In  this  case  the  1  /Q4  behavior 
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Fig.  9.  The  timing  jitter  spectrum  of  an  assumed  spontaneous  emission 
excitation  (solid  line)  compared  to  the  experimental  results  of  [22]  (dashed 
line). 


is  particularly  evident.  All  these  experimental  results  sug¬ 
gest  that  the  source  of  timing  fluctuations  is  not  the  direct 
broadband  excitation  of  timing  jitter  which  would  give  a 
1  /fl2  contribution,  but  rather  a  broadband  source  of  car¬ 
rier  frequency  noise  that  produce  timing  jitter  via  the 
group  velocity  dependence  on  the  carrier  frequency. 

Intensity  noise  and  phase  noise  are  usually  far  from  the 
quantum  limit  in  solid-state  and  gas  lasers.  Unavoidable 
gain  fluctuations  induced  by  fluctuations  in  the  power 
supply  and  length  fluctuations  due  to  environmental  per¬ 
turbations  are  dominant.  On  the  contrary,  as  already 
pointed  out,  when  the  laser  runs  at  the  maximum  of  the 
gain  bandwidth  there  are  no  evident  classical  noise  sources 
for  frequency  fluctuations,  which  are  the  dominant  con¬ 
tribution  to  the  timing  fluctuations  at  low  frequency.  For 
this  reason,  it  is  of  interest  to  evaluate  the  timing  jitter 
spectrum  as  driven  by  quantum  noise.  We  have  gone 
through  a  numerical  example  using  the  following  param¬ 
eters;  2 g  (incremental  gain  per  pass)  =  0.06,  TR  =  1.16 
•  10'8  s,w„  =  2x  1(T7  J,  2D  =  -4000  fs2,  \  =  840 
nm,  Bg  =  1.6  •  1015  s“\  r  =  60/1.763  fs,  noise  en¬ 
hancement  factor  6  =  1  (assuming  that  the  lower  level  of 
the  laser  transition  is  empty).  These  data  have  been  ob¬ 
tained  from  [22]  for  the  Ti:  A1203  laser  described  there. 
The  results  are  shown  in  Fig.  9.  Here,  the  solid  line  is  the 
spectrum  obtained  with  these  data,  the  dashed  line  is  the 
spectral  power  of  the  timing  jitter  obtained  from  [22].  The 
spectra  are  multiplied  by  the  normalization  factor  fl2,  to 
be  consistent  with  the  units  of  dB  over  the  carrier  (dBc) 
usually  reported  in  the  experimental  literature.  The  fre¬ 
quency  dependence  of  the  experimental  results  is  well  re¬ 
produced.  The  peaks  at  50,  150  Hz,  etc.,  which  show  up 
in  the  experimental  spectrum  corresponding  to  amplitude 
noise  in  the  pump  laser,  can  be  attributed  to  low  fre¬ 
quency  fluctuations  of  the  spontaneous  emission.  The  the¬ 
oretical  spectrum  obtained  assuming  only  quantum  fluc¬ 
tuations  as  a  source  of  frequency  noise  is  around  20  dB 
below  the  experimental  spectrum.  This  corresponds  to 
calculated  rms  timing  jitter  one  order  of  magnitude  lower 
than  the  measured  one.  The  discrepancy  may  be  due  to  an 
incomplete  knowledge  of  the  experimental  parameters,  to 
offset  of  the  gain  line  center  from  the  passive  cavity  loss 
spectrum;  to  the  presence  of  Raman  gain  in  the  Ti :  Al203 
which  produces  an  intensity  dependent  shift  of  the  laser 


frequency,  that  couples  intensity  and  frequency  fluctua¬ 
tions;  and/or  to  incomplete  inversion.  Since  the  timing 
jitter  due  to  quantum  noise  is  not  far  from  the  measured 
value,  large  improvement  of  the  timing  jitter  can  be  ob¬ 
tained  only  by  means  of  active  stabilization  of  the  cavity 
[25]  or  by  an  active  “retiming”  of  the  pulse  stream  [22], 
[27  ]. 


IX.  Conclusion 

We  have  presented  an  analysis  of  noise  in  passively 
mode-locked  lasers  with  particular  applicability  to  APM 
or  KLM-locked  lasers  in  which  the  Kerr  nonlinearity 
would  introduce  an  appreciable  chirp  unless  compensated 
by  negative  GVD.  In  this  case,  the  system  obeys  per¬ 
turbed  soliton  equations.  To  the  extent  that  (40)  through 
(43)  are  dynamic  equations  of  motion  with  a  simple  phys¬ 
ical  interpretation,  our  approach  may  predict  noise  behav¬ 
ior  that  is  not  strictly  based  on  soliton  perturbation  theory. 
We  consider  both  classical  and  quantum  noise  sources, 
generally  6  function  correlated  on  the  slow  time  scale  of 
multiple  resonator  transit  times.  When  the  pulsewidth  and 
timing  fluctuations  are  uncorrelated,  a  set  of  symmetric 
spectra)  lines  spaced  by  2t/Tr  is  predicted  in  the  low 
frequency  spectrum  of  the  detector  current.  This  is  in  gen¬ 
eral  agreement  with  von  der  Linde’s  predictions.  In  the 
presence  of  correlation,  these  spectral  lines  may  become 
asymmetric.  The  interference  of  timing  jitter  from  the  low 
frequency  spectra,  possible  for  actively  mode  locked  sys¬ 
tems,  may  lead  to  erroneous  results  when  applied  to  the 
passively  mode-locked  case.  The  spectral  lines  of  the 
photocurrent  show  a  fl~4  dependence  for  fl  >  1  / tp,  where 
7p  is  a  critical  time  related  to  the  gain  bandwidth,  (61). 
These  predictions  are,  generally,  in  agreement  with  ex¬ 
perimental  observations. 
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